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A repeating fast radio burst source in a globular cluster 


F. Kirsten!, B. Marcote?, K. Nimmo**, J. W. T. Hessels*?, M. Bhardwaj>*, S. P. Tendulkar”, A. Keimpema?, J. Yang!, 
M. P. Snelders*, P. Scholz’, A. B. Pearlman$6!®-11.12, C, J, Law!®!4, W. M. Peters!>, M. Giroletti!®, Z. Paragi?, C. Bassa’, 
D. M. Hewitt*, U. Bach!’, V. Bezrukovs!®, M. Burgay!®, S. T. Buttaccio”, J. E. Conway!, A. Corongiu!?, R. Feiler?!, 

O. Forssén!, M. P. Gawroński?! , R. Karuppusamy!’, M. A. Kharinov?”?, M. Lindqvist!, G. Maccaferri!®, A. Melnikov??, 
O. S. Ould-Boukattine*, A. Possenti!®?3, G. Surcis!9, N. Wang”, J. Yuan”, K. Aggarwal*>®, R. Anna-Thomas?>°, 

G. C. Bower”, R. Blaauw’, S. Burke-Spolaor”>”®?8, T. Cassanelli?”’, T. E. Clarke!>, E. Fonseca>®*>?6, B, M. Gaensler??, 
A. Gopinath*, V. M. Kaspi>*°, N. Kassim), T. J. W. Lazio”, C. Leung?!*?, D. Z. Lil’, H. H. Lin?, K. W. Masui?! ?, 

R. Mckinven?, D. Michillić, A. Mikhailov”, C. Ng’, A. Orbidans!®, U. L. Pen??.28:34, E, Petroff+>35, M. Rahman*®, 
S. M. Ransom?’, K. Shin?!**, K. M. Smith**, I. H. Stairs*®, and W. Viemmings! 


‘Department of Space, Earth and Environment, Chalmers University of Technology, Onsala Space Observatory, 439 92, Onsala, Sweden 
Joint Institute for VLBI ERIC, Oude Hoogeveensedijk 4, 7991 PD Dwingeloo, The Netherlands 
3 ASTRON, Netherlands Institute for Radio Astronomy, Oude Hoogeveensedijk 4, 7991 PD Dwingeloo, The Netherlands 
4 Anton Pannekoek Institute for Astronomy, University of Amsterdam, Science Park 904, 1098 XH, Amsterdam, The Netherlands 
> Department of Physics, McGill University, 3600 rue University, Montréal, QC H3A 2T8, Canada 
®McGill Space Institute, McGill University, 3550 rue University, Montréal, QC H3A 2A7, Canada 
7Department of Astronomy and Astrophysics, Tata Institute of Fundamental Research, Mumbai, 400005, India 
8National Centre for Radio Astrophysics, Post Bag 3, Ganeshkhind, Pune, 411007, India 
°Dunlap Institute for Astronomy & Astrophysics, University of Toronto, 50 St. George Street, Toronto, ON M5S 3H4, Canada 
Division of Physics, Mathematics, and Astronomy, California Institute of Technology, Pasadena, CA 91125, USA 
'IMcGill Space Institute (MSI) Fellow 
? ERQNT Postdoctoral Fellow 
Cahill Center for Astronomy and Astrophysics, California Institute of Technology, 1216 E California Boulevard, Pasadena, CA 91125, 
USA 
'4Qwens Valley Radio Observatory, California Institute of Technology, Pasadena, CA 91125, USA 
SRemote Sensing Division, US Naval Research Laboratory, Washington, DC 20375, USA 
‘6Tstituto Nazionale di Astrofisica, Istituto di Radioastronomia, via Gobetti 101, I-40129 Bologna, Italy 
17Max Planck Institute for Radio Astronomy, Auf dem Heugel 69 53121 Bonn, Germany 
'8Engineering Research Institute Ventspils International Radio Astronomy Centre (ERI VIRAC) of Ventspils University of Applied 
Sciences (VUAS), Inzenieru street 101, Ventspils, LV-3601, Latvia 
'Tstituto Nazionale di Astrofisica, Osservatorio Astronomico di Cagliari, Via della Scienza 5, I-09047, Selargius, Italy 
0Tstituto Nazionale di Astrofisica, Istituto di Radioastronomia Radiotelescopio di Noto, C.da Renna Bassa Loc. Case di Mezzo C.P. 161, 
96017, Noto (SR), Italy 
*lTnstitute of Astronomy, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Grudziadzka 5, 87-100 Toru’n, 
Poland 
Institute of Applied Astronomy of the Russian Academy of Sciences, Kutuzova Embankment 10, St. Petersburg, 191187, Russia 
3 Universita di Cagliari, Dipartimento di Fisica, S.P. Monserrato-Sestu Km 0,700 - 09042 Monserrato (CA), Italy 
?% Xinjiang Astronomical Observatory, 150 Science 1-Street, Urumqi, Xinjiang 830011, China 
25 Department of Physics and Astronomy, West Virginia University, P.O. Box 6315, Morgantown, WV 26506, USA 
6Center for Gravitational Waves and Cosmology, West Virginia University, Chestnut Ridge Research Building, Morgantown, WV 26505, 
USA 
27 Academia Sinica Institute of Astronomy and Astrophysics, 645 N. A’ohoku Place, Hilo, HI 96720, USA 
8Canadian Institute for Advanced Research, CIFAR Azrieli Global Scholar, MaRS Centre West Tower, 661 University Ave. Suite 505, 
Toronto ON MSG 1M1, Canada 
David A. Dunlap Department of Astronomy & Astrophysics, University of Toronto, 50 St. George Street, Toronto, ON M5S 3H4, Canada 
3°Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Dr, Pasadena, CA 91109, USA 
3!1MIT Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, 77 Massachusetts Ave, Cambridge, 
MA 02139, USA 
3?Department of Physics, Massachusetts Institute of Technology, 77 Massachusetts Ave, Cambridge, MA 02139, USA 
33Canadian Institute for Theoretical Astrophysics, University of Toronto, 60 St. George Street, Toronto, ON, M5S 3H8, Canada 
34Perimeter Institute for Theoretical Physics, 31 Caroline Street North, Waterloo, ON, N2L 2Y5, Canada 
35Veni Fellow 
36Sidrat Research, PO Box 73527 RPO Wychwood, Toronto, ON M6C 4A7, Canada 
37National Radio Astronomy Observatory, 520 Edgemont Rd., Charlottesville, VA 22903, USA 
8Dept. of Physics and Astronomy, University of British Columbia, 6224 Agricultural Road, Vancouver, BC V6T 1Z1 


Fast radio bursts (FRBs) are exceptionally luminous 
flashes of unknown physical origin, reaching us from 
other galaxies'. Most FRBs have only ever been seen 
once, while others flash repeatedly, though sporadi- 
cally**. Many models invoke magnetically powered 
neutron stars (magnetars) as the engines producing 
FRB emission*>. Recently, CHIME/FRB announced 
the discovery° of the repeating FRB 20200120E, com- 
ing from the direction of the nearby grand design spi- 
ral galaxy M81. Four potential counterparts at other 
observing wavelengths were identified® but no defini- 
tive association with these sources, or M81, could be 
made. Here we report an extremely precise localisa- 
tion of FRB 20200120E, which allows us to associate it 
with a globular cluster (GC) in the M81 galactic sys- 
tem and to place it ~2 pc offset from optical center of 
light of the GC. This confirms® that FRB 20200120E 
is 40 times closer than any other known extragalactic 
FRB. Because such GCs host old stellar populations, 
this association strongly challenges FRB models that 
invoke young magnetars formed in a core-collapse su- 
pernova as powering FRB emission. We propose, in- 
stead, that FRB 20200120E is a highly magnetised neu- 
tron star formed via either accretion-induced collapse 
of a white dwarf or via merger of compact stars in a bi- 
nary system’. Alternative scenarios involving compact 
binary systems, efficiently formed inside globular clus- 
ters, could also be responsible for the observed bursts. 

M81 is one of the most massive nearby galaxies, at 
3.6Mpc®. We targeted the previously best-known posi- 
tion® of FRB 20200120E several times during February 
to May 2021 with an array composed of up to eleven radio 
telescopes that are part of the European Very Long Base- 
line Interferometry (VLBI) Network (EVN; see Methods). 
We observed at a central radio frequency of ~ 1.4 GHz and 
recorded raw voltages from all telescopes. At Effelsberg 
and at the Sardinia Radio Telescope (SRT), we collected 
data in parallel using pulsar data recorders (Methods). 


In total we detected five bursts from FRB 20200120E, 
with dispersion measures close to the previously reported ° 
DM = 87.8pccm™*. Two bursts were detected on 2021 
February 20 (called B1 and B2 below), two bursts on 2021 
March 7 (B3 and B4), and one burst on 2021 April 28 
(B5). Bursts B2-B5 were found by blindly searching both 
the Effelsberg voltages and PSRIX° data (Methods), while 
B1 was only detected in the voltage data as it occurred 
outside the recording times of the PSRIX instrument !?. In 
Figure 1 we show the dedispersed dynamic spectra and 
frequency averaged burst profiles. Burst fluences range 
from 0.13 —0.71 Jy ms and total burst envelopes span only 
~ 100 — 300 us (Table 1). A detailed, ultra-high-time- 
resolution analysis of the burst properties and energetics 
is presented in a companion paper!”, 


Correlation of the data, in order to produce ‘visibilities’ 
for interferometric imaging, was performed at the Joint In- 
stitute for VLBI ERIC. To achieve the best-possible sensi- 
tivity, we used the coherent-dedispersion mode of the soft- 
ware correlator SFXC!!, applying DM = 87.77 pccm™, 
which we derived from a manual inspection of the bursts 
(Methods). 

After an initial rough localisation via delay mapping 
(Methods, accurate to several arcseconds), we individu- 
ally imaged the five bursts, where each data set spans 
only the width of each burst in time (Table 1). Given 
the snapshot nature of the correlations, the rather sparse 
arrays in each run, and the fact that the bursts only cov- 
ered a fraction of the observed bandwidth, the images 
from the individual bursts result in elongated fringe pat- 
terns, hindering an individual, unambiguous localisation 
of each burst at the level of the synthesised beam size 
(see Figure 2a-d). Therefore, we created a data set that 
is the combination of the visibilities from all bursts ex- 
cept BI — as it was too faint to produce a useful im- 
age, and we therefore exclude it from the localization 
analysis. These data allowed us to unambiguously pin- 
point the position of FRB 20200120E in the field (see Fig- 
ure 2e,f). The derived coordinates of FRB 20200120E 
in the International Celestial Reference Frame (ICRF) !” 
are RA (J2000) = 9557™54.69935° + 1.2 mas Dec (J2000) 
= 68°49’0.8529” + 1.3 mas (see Methods). These coor- 
dinates coincide with the location of the globular cluster 
[PR95] 30244 '?, which is part of the M81 globular cluster 
system |*, 

Figure 3a shows the position of FRB 20200120E with 
respect to [PR95] 30244 in a combined three-colour Sub- 
aru image made with i’, r’, and g’ filters mapped to red, 
green, and blue channels, respectively. The galaxy at 
the bottom-left of [PR95] 30244 is a background SDSS 
galaxy at redshift z = 0.194, i.e. at a z much larger 
than the maximum possible® value (z < 0.03) based on 
FRB 20200120E’s dispersion measure. We performed ra- 
dial fits to the brightness distribution of [PR95] 30244 us- 
ing a Moffat profile! in both RA and Dec for all three 
bands after fitting and subtracting a bilinearly varying 
background (to account for the presence of the back- 
ground galaxy). The average position of the centroid of 
[PR95] 30244 is RA (J2000) = 9"57"54.7135° + 7 mas, 
Dec (J2000) = 68°49’0.766” + 4 mas (statistical), well 
in agreement with previous measurements !? and the posi- 
tion of the source in the Gaia Early Data Release 3 Cata- 
logue !6!7 (Methods), in which positions are well aligned 
with the ICRF”. 

The centre of the FRB localisation is ~ 116mas off- 
set from the optical centre of light of [PR95] 30244 (Fig- 
ure 3b, corresponding to ~ 54% of its effective radius, 
Methods). Given the astrometric uncertainty of the FRB 


localisation (~ 1.3 mas) and the optical image registra- 
tion error with respect to the ICRF (~ 15 mas; see Meth- 
ods), we conclude that FRB 20200120E is located sig- 
nificantly (> 7o confidence level) offset from the opti- 
cal centre of light of [PR95] 30244. The optical angu- 
lar size of [PR95] 30244 (0777; Methods) in combina- 
tion with FRB 20200120E’s offset from M81 (19.6’) and 
the number of GCs predicted '® to be part of the galaxy 
(300 + 100), allow us to estimate the probability of chance 
alignment Pee < 1.7 x 1074 (Methods). From such a 
very low Pec value, we conclude that the association of 
FRB 20200120E and [PR95] 30244 is robust. 

Figure 4a shows a deep continuum map that was cre- 
ated from the combination of the data of all VLBI obser- 
vations. We find no persistent source at the location of 
FRB 20200120E above a 5c confidence level (rms back- 
ground noise level of 10 wy beam~'). Also shown in Fig- 
ure 4 are continuum images obtained with the Karl G. Jan- 
sky Very Large Array (VLA) at 1.5 GHz and at 340 MHz 
between December 2020 to January 2021 (Methods). 
These maps have noise levels of 6.5 wy beam™! and 
320 uJy beam™!, respectively. Also here, no persistent 
source is detected at the position of FRB 20200120E in 
either of the images. For a 1.5GHz radio flux density 
limit of 20 wy (30) and a distance of 3.6 Mpc, we limit 
the radio luminosity L, < 3.1 x 10? ergs"! Hz7!. This lu- 
minosity limit is ~ 10° times lower than that of any other 
extragalactic FRB '° and almost 10° times lower than the 
radio luminosity of the persistent source in the vicinity of 
FRB 20121102A”. 

We find no evidence of an X-ray source at the loca- 
tion of FRB 20200120E in archival Chandra observations 
(Methods). This results in a 0.5—10 keV luminosity upper 
limit of 2 x 10°’ ergs"! (37) at the distance of 3.6 Mpc. 
A detailed analysis of ongoing X-ray follow-up obser- 
vations of the region will be presented in Pearlman et 
al. (in prep.). Similarly, no sources are reported at the 
location of M81 in any of the Fermi-LAT catalogues”!. 
The nearest catalogued source (at a separation of 52’) is 
4FGL J0955.7+6940, known to be associated with M82”. 

Within the context of FRB models that invoke a young, 
highly magnetized neutron star (NS) powered primarily 
by the decay of its magnetic field*”, i.e. a magnetar, it 
is hard to reconcile FRB 20200120E’s association with an 
old globular cluster using the standard core-collapse su- 
pernova formation channel of magnetars. Instead, because 
of their extreme stellar densities, globular clusters are 
known to form short-orbital-period binaries at a high spe- 
cific rate?>~*>, We thus propose that FRB 20200120E is a 
magnetar formed via accretion-induced collapse (AIC)*° 
of a white dwarf (WD) or via merger-induced collapse 
(MIC) of a WD-WD, NS-WD or NS-NS binary”? — 
systems that are common in globular clusters and, like 


FRB 20200120E, are found concentrated towards their 
core3? (Methods). The lack of a persistent radio or X-ray 
source at the position of FRB 20200120E is expected in 
an AIC/MIC scenario, as any emission generated during 
collapse fades on short time scales (< 1 yr)’. 

The globular cluster host of FRB 20200120E also sug- 
gests some alternatives to the magnetar class of FRB mod- 
els. FRB 20200120E could be a compact binary system 
— such as a tight WD-NS system in a pre-merger phase 
or a magnetised NS with a planetary companion*!*? — 
in which the bodies are interacting magnetically. Simi- 
larly, a binary millisecond pulsar with a strong magnetic 
field formed via AIC and that was subsequently spun-up 
via accretion**** could act as an FRB engine. Such a sys- 
tem could also be observable as a low mass X-ray binary 
(LMXB)*>, as would an accreting black hole (BH). In 
such an LMXB model, the radio bursts could be generated 
via magnetic reconnection in a relativistic jet or where the 
jet shocks with the surrounding medium and creates a syn- 
chrotron maser*°. Except for the most luminous LMXBs 
(Lx = 10% erg s71), our observations cannot rule out such 
systems. However, none of the ~ 200 Galactic LMXBs 
has been seen to generate FRBs. In some cases, ultra- 
luminous X-ray sources (ULXs) 37 have been shown to be 
NSs accreting at hyper-Eddington rates**, though some 
may be systems with a more massive BH primary”. We 
note that ULXs have been associated with extragalactic 
globular clusters*” but such systems are ruled out by our 
X-ray limit unless their luminosity varies in time by more 
than two orders of magnitude. Additionally, the associa- 
tion with a globular cluster rules out a high mass X-ray bi- 
nary origin of FRB 20200120E and the projected offset of 
x 2 pc from the centre of light of [PR95] 30244 excludes 
the association of FRB 20200120E with, e.g., a massive 
(= 10M.) stellar mass BH or a putative intermediate mass 
black hole at the core of [PR95] 30244. 

The association of FRB 20200120E with a globular 
cluster adds to the diversity of environments in which re- 
peating FRBs have been found. While FRB 20121102A 
resides in a dwarf galaxy*!, the host of FRB 20180916B 
is a spiral galaxy!’ and FRB 20201124A was localised 
to a massive star forming galaxy**. Previously localised 
repeaters have been associated with nearby star forming 
regions **~*°, favouring the core-collapse supernova chan- 
nel for the formation of a young magnetars, as the rate of 
AIC and MIC is much lower. The lack of a persistent radio 
source for all but FRB 20121102A may suggest a range in 
the possible ages of such magnetars. In a globular cluster 
environment, however, the recent core-collapse of a mas- 
sive star is very unlikely. Thus, this suggests a diversity in 
formation channels for magnetars as FRB engines. 
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Figure 1: Dispersion-corrected time series and dynamic spectra of the five FRB 20200120E bursts Frequency- 
averaged time series of the bursts are displayed in panels a-e. The top right corner of each panel shows the burst 
name used in this work. The blue-coloured bars highlight the +20 burst width used to measure the burst fluence. The 
dynamic spectra of each burst are shown in panels f-j. The red marks represent the edges of the subbands. Data that 
have been removed due to contamination by radio frequency interference have not been plotted. In all panels the data 


are plotted with 32 us and 2 MHz time and frequency resolution, respectively (with the exception of B1 and B5 which 
are plotted with 4 MHz frequency resolution). 


Table 1: Burst properties. 


Fluence®* Peak Flux Width? Gate width ° 
Burst MJD ° [Jy ms] Peak S/N‘ Density>* [Jy] [us] [us] 
Bl 59265.88304437179 0.13 0.03 7.8 0.9 + 0.2 1561 290 
B2 59265.88600912486 0.63 0.12 54.9 6.6+ 1.3 62+ 1,93+0.5! 150 
B3 59280.69618745651 0.52 +0.10 64.5 78+ 1.6 46.7+0.1 126 
B4 59280.80173397988 0.71 0.14 47.0 5.7+41.2 117+1 386 
B5 59332.50446581106 0.09 + 0.02 11.6 1.4+0.3 56.6 + 0.1 173 


a Corrected to the Solar System Barycentre and to infinite frequency assuming a dispersion measure of 87.75 pe cm 
reference frequency 1502 MHz and dispersion constant of 1/(2.41x 107+) MHz? pc~! cm? s. 
The times quoted are dynamical times (TDB). 
> The receiver temperature of Effelsberg is 20 K and the telescope gain is 1.54 K Jy~!. We additionally consider 
a sky background temperature of 0.8 K, by extrapolating from the 408 MHz map“, 
using a spectral index of —2.7 8, and 3 K from the cosmic microwave background 
49. We take a conservative 20 % error on these measurements, arising due to the uncertainty 
in the Effelsberg receiver temperature and gain. 
e Computed using the frequency range over which the burst is bright. 
“Defined as 1/ v2 multiplied by the full-width at half-maximum of the autocorrelation function (ACF). 
Note we use the +20 burst width region to determine the burst fluence (see Nimmo et al. (submitted) for details). 
e Width of the gate used for the interferometric correlation of each burst. For B2 the gate was centered on the first, brighter 
component to maximise S/N. 
Width per burst component. 
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Figure 2: Localisation plots for FRB 20200120E. Normalised dirty images of the individual bursts (B2, B3, B4 and 
B5; a-d) along with a dirty image of the four bursts combined (e), produced by applying a natural weighting to the 
data. For visualisation purposes we clip the colour scale at zero, i.e. only positive peaks are displayed. The white 
circles are centred on the location of the globular cluster [PR95] 30244 as derived from the Subaru image in Figure 3. 
Their size indicates the region that contains 90% of the globular cluster’s optical emission. The synthesized beams of 
each image are displayed as grey ellipses in the bottom left corner of each panel. f: The cleaned image and resulting 
localisation of FRB 20200120E, as derived from the combined data sets of four bursts. The resulting coordinates 
of FRB 20200120E (highlighted by the white marker) are RA (J2000) = 9°57™54.69935° + 1.2 mas, Dec (J2000) 
= 68°49’0.8529” + 1.3 mas. 
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Figure 3: Optical images of the FRB 20200120E host and surrounding field. Left (a): 40”x 40”g', r', and i’ 
band image of [PR95] 30244 acquired with Hyper Suprime-cam. The small red ellipse is centred at the location of 
FRB 20200120E. In panel b we show the zoomed-in r’ band image of [PR95] 30244. The grey circle represents the 
estimated position of the centre of [PR95] 30244 and its 100 uncertainty (dominated by the optical-to-radio reference 
frame tying). The small red ellipse is the same as in panel a, and also represents the 100 positional uncertainty region 
of FRB 20200120E. Panels c and d show cross-sections of the brightness distribution of the cluster (blue solid lines) 
with the Moffat profile that we fit overlaid in black. Indicated in solid grey lines are the PSFs as measured from stars 
in the images. Note that scatter in the PSFs is smaller than the linewidth. 
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Figure 4: Continuum maps of the field around FRB 20200120E. a: 1.4-GHz EVN continuum image after com- 
bining the three epochs (EK048B, EK048C, and EK048F); b: 1.5-GHz Realfast, and c: 340-MHz VLITE continuum 
image. The red circles indicate the 100 (for EVN) and 10000 (for Realfast, VLITE) positional uncertainty region of 
FRB 20200120E. Note the very different scales between the three panels. We clip all values below zero and above 
60 wy beam™! (EVN), 50 uJy beam™! (Realfast) and 3 mJy beam™! (VLITE) for visualisation purposes. The black 
ellipse in the bottom left corner of each panel indicates the synthesised beam size and position angle. 


Methods 


Observations and Data reduction 
VLBI observations 


As part of our ongoing very long baseline interferome- 
try (VLBI) campaign called PRECISE (Pinpointing RE- 
peating ChIme Sources with the EVN), we observed the 
field centred on the coordinates from the discovery® of 
FRB 20200120E: RA (J2000) = 09°57™56.688°, Dec 
(J2000) = 68°49’31.8”. The quoted 90% confidence inter- 
valí of radius ~ 1.5’ is well covered by the field of view 
of the ad-hoc interferometric array we used to observe, 
which is comprised of dishes that are also part of the EVN. 
We observed the field 15 times between February and May 
2021 (Table 2), with burst detections in the data that were 
taken on 2021 February 20 UT 17:00-22:00 (project code 
PR141A), 2021 March 7 UT 15:45-20:45 (PR143A), and 
2021 April 28 UT 11:00-22:00 (PR158A). We observed in 
the 21-cm band (~ 1.4GHz) with slightly different array 
setups in each run. Depending on the capabilities at each 
station, we recorded either 128 or 256 MHz of total band- 
width, divided into 8 or 16 subbands of 16 MHz each. The 
participating stations and their respective frequency cov- 
erage are summarised in Table 3. We ran regular phase- 
referencing observations with a cycle time of 7.5 min; that 
is, 5.5 minutes on target and 2.0 minutes on the phase cal- 
ibrator (J0955+6903, at ~ 0.32° separation). In total, we 
spent 2.93 hrs on the target field in each of PRI41A and 
PR143A, and 6.73 hrs in PR158A. The calibrator source 
J1048+7143 served as fringe finder and bandpass calibra- 
tor and was observed twice for 5 minutes in each run. For 
verification and single-dish time-domain calibration pur- 
poses we observed the pulsar B0355+54 for 2 minutes in 
PR141A and PR143A. In PRI58A the pulsar B2255+58 
was observed for 5 minutes for the same reason. 

We recorded raw voltages (‘baseband’ data, dual circu- 
lar polarisation, 2-bit sampling) in VDIF® or Mark5B>! 
format at all stations. At Effelsberg we also recorded total 
intensity filterbanks with the PSRIX pulsar backend? in 
parallel to the voltage data. These data span the frequency 
range of 1255-1505 MHz with a time and frequency reso- 
lution of 102.4 us and 0.49 MHz, respectively. Similarly, 
at SRT we also recorded in parallel baseband data and to- 
tal intensity filterbanks (with the local Digital Filter Bank 
system, DFB). These DFB-filterbanks have a time and fre- 
quency resolution of 128 us and 1 MHz, respectively, cov- 
ering the frequency range of 1140.5—2163.5 MHz out of 
which we search the usable range of 1210.5—1739.5 MHz. 

The data from both Effelsberg and SRT were searched 
for bursts in two independent pipelines. The baseband 
data were processed with the pipeline outlined in Kirsten 
et al.°2, which converts the raw voltages to filterbank 


format (in this case, total intensity with time resolution 
64 us and frequency resolution 125kHz) and searches 
those with Heimdal1 within +50 pc cm”? of the expected 
DM = 87.818 pe cm™?, as found by Bhardwaj et al.°. The 
burst candidates are classified as either radio frequency in- 
terference (RFI) or potential FRBs using the neural net- 
work classifier FETCH**. The filterbanks as recorded 
with the respective pulsar backends were searched with 
a pipeline that uses the PRESTO suite of tools**°° and a 
classifier based on work by Michilli et al. 56. 


Correlation of the data was performed at the Joint In- 
stitute for VLBI ERIC (JIVE; in the Netherlands) un- 
der proposal EK048 (with codes EK048B for PRI41A, 
EK048C for PR143A, and EK048F for PR158A). In total, 
we ran three correlator passes using the software correla- 
tor SFXC!!: in the first run we correlate all scans con- 
taining the calibrator sources with a standard 2-s time in- 
tegration and 64 channels per 16-MHz subband. A sec- 
ond correlation pass was performed only on the data con- 
taining bursts (their arrival times were determined via the 
search described above), with a higher frequency resolu- 
tion (8 192 channels per subband) and time resolution (the 
correlation gates around the bursts were chosen manually 
to optimise the signal-to-noise for each burst, Table 1). 
The strongest bursts allowed a direct fringe fit on their 
data, which provided a first estimate of their sky position 
to a level of a few arcseconds: RA (J2000) ~ 09557™54.8°, 
Dec (J2000) ~ 68°49’03”. We then used this position to 
re-correlate the burst data at the same frequency resolution 
as in the first correlation pass. This allowed us to directly 
apply the calibration performed on the first pass (contain- 
ing the calibrator sources) to the data containing bursts. 


The full data reduction was performed using the Astro- 
nomical Image Processing System AIPS, >’ and Dirmap*® 
following standard procedures. A-priori amplitude cali- 
bration was performed using the known gain curves and 
system temperature measurements recorded at each sta- 
tion during the observation. These data were not avail- 
able for Torun and Irbene in EK048B, and for Urumqi in 
EK048C. We used nominal system equivalent flux density 
values and flat gain curves to perform the amplitude cali- 
bration in these cases. In EKO48C, the data from Urumqi 
were flagged due to the impossibility to properly convert 
the locally-recorded linear polarizations into circular ones 
(as in the rest of the antennas). Ionospheric dispersive 
delays were calculated from maps of total electron con- 
tent provided by the global positioning system satellites 
and removed via the TECOR task in AIPS. We first cor- 
rected for the instrumental delays and bandpass calibra- 
tion using the calibrator source J1048+7143, and there- 
after fringe-fit the data using all calibrator sources. The 
phase calibrator (J0955+6903; displaying a peak bright- 
ness of ~ 73mJybeam™!) was then imaged and self- 


calibrated to improve the final calibration of the data (for 
both phases and amplitudes). A common source model for 
the phase calibrator was used to improve the calibration of 
both epochs. The obtained solutions were then transferred 
to the burst data, which were finally imaged. 

In the third, and final, correlator pass we processed the 
entire observation, correlating all scans from the target 
with the same time and frequency setup as in the first cor- 
relator pass, in order to obtain a deep continuum image of 
the field. 


VLA 


We observed the field of FRB 20200120E with the Karl G. 
Jansky Very Large Array (VLA) as part of a program to 
localise repeating FRBs identified by CHIME/FRB (pro- 
gram 20B-280). The field was observed in five 1-hour 
blocks from 2020 December 29 through 2021 January 
18. Data were recorded in parallel at both 1.5GHz and 
340 MHz. The VLA antennas were arranged in the A con- 
figuration, which provides baseline lengths up to 30km 
and a typical spatial resolution of 1.3” and 6.0” at 1.5 GHz 
and 340 MHz, respectively. The total on-target integration 
time amounts to 200 minutes. 


Realfast Visibility data from the VLA observations 
were recorded with a sampling time of 3 seconds from 
1 to 2 GHz. Simultaneously, a copy of the data with sam- 
pling time of 10 ms was streamed into the Realfast sys- 
tem*’. We used Realfast to search this data stream for 
FRBs in real time. The typical 1o sensitivity of the VLA 
is 5 mJy beam“! in 10 ms. 

The standard (slow) visibility data were analyzed to 
search for persistent emission associated with the FRB lo- 
cation. We calibrated the data with the VLA calibration 
pipeline (version 2020.1) using 3C147 as a flux calibrator. 
Calibrated visibilities for all five epochs were combined 
and imaged with CASA (version 6.1). 

We imaged the data with robust weighting of 0.5 and 
removed baselines shorter than 1002 to reduce the effects 
of RFI. This produced an image with a synthesized beam 
size of 2” by 1” with position angle of 126°. The sensi- 
tivity in the combined image is 6.5 wy beam™!, which is 
consistent with expectations given the usable bandwidth 
of 400 MHz. 


VLITE The VLA Low-band Ionosphere and Transient 
Experiment VLITE; 60.61 ig a commensal instrument on 
the VLA that records and correlates data from a 64-MHz 
subband at a central frequency of ~ 340 MHz. It operates 
on up to 18 antennas during nearly all regular VLA opera- 
tions. All VLITE data were processed within the VLITE- 
Fast GPU-based real-time system to search the incoming 


voltage stream for dispersed transients®’. For imaging 
purposes, primary calibration and editing for each day of 
visibility data were carried out with the automated VLITE 
processing pipeline. Due to radio interference from satel- 
lites at the upper end of the band, the final usable band- 
width was 38.2 MHz centred at 340.85 MHz. The data 
were subsequently combined, imaged and self-calibrated 
in amplitude and phase using the Obit task MFImage®. 
In order to reduce artefacts from the bright extended radio 
galaxy M82 located ~ 1° northwest of the target position, 
baselines shorter than 4.0kA were removed at this point. 
The final image was created in WSClean™, and corrected 
for the offset primary beam response of VLITE. The im- 
age has an rms of 320 wy beam™!, and a beam of 10.1” 
by 3.6” at a position angle of 132°. 


Archival optical and high energy data 


Hyper Suprime-Cam The field around M81 has been 
well-observed over the years by multiple telescopes. We 
retrieved archival data from the Hyper Suprime-Cam 
on the 8.1-m Subaru telescope® using the SMOKA 
interface. We chose images with seeing better than 
077. We processed the g’, r’, and i’ band images with 
the hscpipev8.4 pipeline®. The pipeline uses the 
PanSTARRs catalogue PS1; ® as an astrometric and pho- 
tometric reference. The typical astrometric residuals were 
~ 50-60 mas, which is the uncertainty we have assumed 
to tie the optical and the radio reference frames. 


Gaia [PR95] 30244 also appears in the Gaia 
Early Data Release 3 Catalogue'®!’ with Source 
ID = 1070264274879949 184, and position RA 
(J2000) = 9°57™54.71402° + 1.6mas, Dec (J2000) 
= 68°49’0.7775” + 1.7 mas. This position is consistent 
(within < 30 confidence level) with the one we have 
derived from the Hyper Suprime-Cam data. The observed 
offset allowed us to estimate the possible systematic un- 
certainties in the optical image registration error (15 mas; 
by adding in quadrature the observed offset between the 
the Gaia position and the one we determined, plus the 


uncertainties on the positions). 


Chandra X-ray Observatory Several deep archival X- 
ray observations are available for the field around M81 
from XMM and Chandra. We selected the archival ob- 
servation with the longest exposure time that covers the 
location of FRB 20200120E, a 26 ks Chandra observation 
Obs. ID 9540, taken with ACIS in FAINT mode, © to 
probe for an X-ray source. The data were reduced us- 
ing CIAO version 4.12 following standard procedures. 
As the source was located about 14’ off-axis, events were 


Table 2: Time ranges of PRECISE runs targeting FRB 20200120E between February and May 2021 


Observation EVN 

Project code Project code Start MJD Stop MJD 
PR141A? EK048B 59265.708 59265.916 
PR143A? EK048C 59280.656 59280.864 
PRI44A 59283.792 59284.000 
PRI45A 59289.750 59289.958 
PRI46A 59295.667 59295.875 
PRI53A 59314.887 59314.972 
PRI58A? EK048F 59332.458 59332.916 
PRIS9A 59336.708 59337.000 
PRI60A 59341.833 59342.072 
PRI61A 59344.771 59344.875 
PR162A 59346.646 59346.895 
PRI63A 59347.417 59347.625 
PR164A 59351.917 59352.166 
PR165A 59358.917 59359.166 
PR166A 59360.708 59360.916 


Epoch with detection 


Table 3: Setups at the different stations during observations used in the analysis 


Telescope Frequency coverage [MHz] Station project code EVN project code 
Effelsberg (Ef) 1254 — 1510 94-20 EK048B/C/F° 
Medicina (Mc) 1350 — 1478 44-20 EK048B/C/F* 
Noto (Nt) 1318 — 1574 44-20 EK048 C 
Irbene (Ir) 1382 — 1510 — EK048B/C/F 
Torun (Tr) 1254 — 1510 DDT* EK048B/C/F° 
Westerbork (Wb) 1382 — 1510 DDT? EK048B/C/F 
Urumqi (Ur) 1382 — 1510 DDT? EK048 C/F 
Sardinia (Sr) 1360 — 1488 44-20 EK048 F 
Onsala (O8) 1382 — 1510 DDT? EK048B F’ 
Badary (Bd) 1382 — 1510 DDT? EK048 F 
Svetloe (Sv) 1382 — 1510 DDT? EK048 F 
Zelenchukskaya (Zc) 1382 — 1510 DDT? EK048 F 
VLA-VLITE 320 — 384 20B-280 — 
VLA-Realfast 1300 — 1700 20B-280 — 


*Director’s Discretionary Time 


Only stations recording the part of the band where the burst (B5) was detected with significant emission. 
“Only one subband overlapping with the part of the band where the burst (B5) was detected with significant emission. 


extracted in a large 10’-radius region centred on the po- 
sition of FRB 20200120E. We also extract events from a 
60’’-radius region away from the source at a similar off- 
axis angle to estimate the background count rate. The X- 
ray count-rate in the source extraction region, 4.4 x 1076 
counts s~! arcsec~2, is consistent with the background re- 
gion rate of 3.9x 1076 counts s~! arcsec”. To place a limit 
on an X-ray source at the location of FRB 20200120E we 


use the Bayesian method of Kraft et al.’°, which results 
in a 0.5-10 keV source count rate upper limit of 1 x 1077 
counts s~! (37). Taking into account the spectral response 
for the off-axis location of the source (via an ancillary re- 
sponse file created by the CIAO tool specextract), and 
assuming a photoelectrically absorbed power-law source 
spectrum with a spectral index of I’ = 2 and a hydrogen 
column density of Ny = 102! cm~”, this count rate limit 


corresponds to a 0.5—10 keV absorbed flux upper limit of 
1 x 107'4ergcem™?s7!. 


Fermi-LAT The Large Area Telescope (LAT) onboard 
the Fermi satellite provides a uniform sensitivity survey 
of the whole sky in the energy range between 100 MeV 
and 300 GeV. We searched all of the publicly available 
catalogues for counterparts up to the latest published”! re- 
lease, 4FGL-DR2, with null results. However, even the 
most luminous known Galactic globular cluster (Terzan5), 
whose luminosity is (42.4 + 1.5) x 10*4 erg s7! in the 0.1- 
100 GeV energy range7?, would have a gamma-ray flux 
of only (2.00 + 0.07) x 107'® erg cm~? s7! at the distance 
of M81. This is nearly three orders of magnitude dimmer 
than the faintest source detected in the 4FGL-DR2 cata- 
logue. 


Analysis 
Dispersion measure refinement 


To refine the DM for further analysis, we maximised the 
signal to noise ratio (S/N) of a very narrow spike in B3!° 
(Figure 1g) and find DM = 87.7527 + 0.0003 pe cm™?. 
This value is comparable to, but formally deviates from 
DM = 87.818 + 0.007 pe cm? found by Bhardwaj et al. ê, 
where a weighted average of three bursts was used. Possi- 
ble explanations for this discrepancy are the lack of short- 
timescale structure in the CHIME/FRB bursts°, the po- 
tential for non-dispersive time-frequency drifting”? or a 
time-variable DM. 


Milliarcsecond localization of FRB 20200120E 


We imaged both the individual bursts separately, as well 
as a data set produced by the combination of all individ- 
ual burst visibilities. Figure 2a-d display the dirty maps 
(i.e. the inverse Fourier Transform of the visibilities with- 
out applying any deconvolution, i.e. “CLEANing’) of the 
bursts that were detectable in the correlated data, using a 
natural weighting of the data. B1 was too faint to produce 
a useful image, and we therefore exclude it from the lo- 
calisation analysis. B5 was only detected in the lower half 
of the observed band, where most of the antennas were 
not recording (see Table 3). We therefore only used data 
from this part of the band. Figure 2e shows the dirty map 
of the combined data of the visibilities from all bursts. 
In this map we obtained an emission pattern consistent 
with the one expected from the dirty beam (the inverse 
Fourier Transform of a point-like source), allowing us to 
unambiguously identify the position of FRB 20200120E. 
The observed emission reached a 120 confidence level, 
with the secondary sidelobes in the fringe pattern being 


~ 66% of the peak emission. This provided a robust lo- 
calization in the map, as it would require a noise fluctu- 
ation of = 70 to produce such peak emission. We also 
conducted different approaches during the imaging of the 
data: different weighting schemes, and selecting different 
subsets of antennas in the array. The derived position of 
FRB 20200120E was robust across all these approaches. 
Figure 2f displays the final, ‘CLEANed’, image of the 
combined bursts. 

The final coordinates of FRB 20200120E are RA 
(J2000) = 9°57™54.69935° + 1.2mas, Dec (J2000) = 
68°49'0.8529” + 1.3 mas. We note that the quoted uncer- 
tainties reflect the statistical uncertainties from the mea- 
sured position of FRB 20200120E (0.7 and 0.4mas in 
RA and DEC, respectively), the uncertainties in the ab- 
solute International Celestial Reference Frame position of 
the phase calibrator (J0955+6903; 0.11 mas), and the sys- 
tematic uncertainty associated with the phase-referencing 
technique ’* of ~ 0.9 and 1.2 mas, in RA and DEC, respec- 
tively. 

We combined the continuum VLBI data from the three 
epochs to produce a deep image of the field around 
FRB 20200120E to search for persistent emission. No sig- 
nificant sources above a 5o confidence level (with an rms 
of 10 uJybeam™!) are detected on milliarcsecond scales 
(Figure 4a). In the VLA data taken at 1.5 GHz (Fig- 
ure 4b), we did identify two sources with peak bright- 
ness of 110 wy beam™! and 73 wy beam, and offset by 
6” and 9”, respectively. The projected density of ra- 
dio sources of this brightness is roughly 1000 to 3000 
per square degree’. Therefore, we expect between 1- 
3 sources within 1’ of FRB 20200120E by chance. The 
closer of the two nearby radio sources is within 0.2” of 
a PS1 source with i = 21.3 mag. This host galaxy has 
a photometric redshift in the PS1 STRM catalogue’® of 
0.67 + 0.2, so it is most likely a background galaxy. The 
other identified radio source has no PS1 counterpart. We 
note that neither of these two sources exhibit significant 
compact emission on milliarcsecond scales. 


[PR95] 30244 and chance coincidence probability 


We measured the full-width at half maximum (FWHM; 
“seeing”) of the coadded i’, r’, and g’ images of 
[PR95] 30244 to be 063, 0/57, and 0/62, respectively, 
using profile fits to bright, isolated stars. The bright- 
ness distribution of [PR95] 30244 has an FWHM of 077, 
0’70, and 075 in the same three images. Subtracting the 
FWHMs of the isolated stars in quadrature, we estimate 
that the intrinsic FWHM of [PR95] 30244 is about 0742, 
corresponding to about 7.4 pc at a distance of 3.63 Mpc. 
To estimate the probability of chance coincidence” for 
an M81 globular cluster in the FRB localisation region, 


we use a circular localisation region with a radius (R) = 
max{2Reg of [PR95] 30244, maximum seeing-limited size 
of [PR95] 30244 in the Hyper Suprime-images} = 0°77. 
Perelmuter & Racine '* parameterized the projected areal 
number density (0Gc/ arcmin”) of the M81 globular clus- 
ters as a function of their angular offset from M81 (r, in 
arcmin) as, log )(eGc) = —2.07 xlogj9(r) + 0.82+0.05, for 
0 < log, (r) < 1.4. At the offset of FRB 20200120E from 
M81 (19.6’), pac = 0.014 arcmin™? x 3.8x10~° arcsec”. 
Assuming a Poisson distribution of M81 globular clus- 
ters at r = 19.6’, the probability of finding at least one 
globular cluster by chance within a radius R is given by 
Poo = 1- exp(—71R*pgc) x 7x 1076. A more conser- 
vative estimate of Pec can be derived by assuming that 
all the predicted 300 + 100 M81 globular clusters !8 are 
uniformly distributed within the angular area of radius = 
19.6’ so that pgc = 4.6-9.2 x 1075, and consequently 
Peo = 0.85-1.7 x 1074. From such a very low Pec value, 
we conclude that the association of FRB 20200120E and 
[PR95] 30244 is robust. 


Modelling of [PR95] 30244 


To estimate important physical properties of 
[PR95] 30244, such as its stellar mass, metallicity, 
stellar population age and V-band extinction, we used the 
Prospector code? for stellar population inference. 
We modelled the SDSS photometry (Table 5, from the 
SDSS DR12 catalogue *”) of [PR95] 30244 and fit a 
five-parameter (Table 6) ‘delayed-tau’ model for the star 
formation rate SFR(t) œ texp(—t/T), where t is the time 
since the formation epoch of the galaxy, and t is the 
characteristic decay time of the star-formation history 
of [PR95] 302448!8?, As we are modelling a globular 
cluster, we did not include nebular line emission but 
only enabled a dust emission model? in our fitting. 
The best-fit spectral energy distribution (SED) profile of 
[PR95] 30244 is shown in Figure 6. Prospector also 
enables Markov Chain Monte Carlo (MCMC) sampling 
of the posterior to estimate uncertainty in the best-fit 
values of the physical properties of [PR95] 30244. We 
show the corner plot of the MCMC-analysis in Figure 7 
and list the results in Table 4. Using the relation 84, 
[Fe/H] = 1.024 x log(Z/Zo), we estimate the [Fe/H] of 
[PR95] 30244 = — 1.837086, which is in good agreement 
with earlier estimates !**°. We can get an estimate of the 
velocity dispersion (o+) of [PR95] 30244 using the virial 
theorem: or ~ V2GM/3Re, where G is the gravitational 
constant. Using the M and Rep values of [PR95] 30244 
from Table 4, we estimate o, ~ 22kms7!. 
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Possible MIC-models 
FRB 20200120E 


for the formation of 


The most likely of the MIC-models is that of a merging 
WD-WD system as those dominate the cores of globu- 
lar clusters®°, while a NS-NS progenitor system is less 
probable*’. A typical globular cluster with a total mass 
of ~ 2 x 10° Mo can host 10 — 20 MSPs formed via AIC 
or MIC**. At birth, such NSs would be extreme objects 
with high rotation rates and strong magnetic fields, po- 
tentially capable or generating FRBs. We note that the 
metallicity limit (log[Fe/H] > —0.6) set for the forma- 
tion of young NSs in globular clusters®* is much higher 
than the metallicity that we estimate for [PR95] 30244 
(log[Fe/H] = —1.83, Table 4). Nevertheless, non-recycled 
pulsars have been found towards some globular clus- 
ters**-°°. These young pulsars are formed at a compara- 
ble rate to MSPs in globular clusters, but they have much 
shorter active lifetimes (Myr versus Gyr). 


Constraints on the Galactic halo DM and RM 


Since FRB 20200120E is extragalactic, its disper- 
sive delay is probing the full extent of the Milky 
Way ISM (DMX) and halo (DMM) along this 
line of sight. There are also extragalactic contribu- 
tions from the intergalactic medium DMjgm and the 
host galaxy (DMmg1), leading to the total observed 
DM = DMN + DMhy + DMicm +DMmsi. The inter- 
galactic medium between the Milky Way and M81 con- 
tributes on the order’! of DMigm ~ 1 pe cm. To es- 
timate DMygi, we only consider the contribution from 
the halo of M81, as both the M81 disk and the glob- 
ular cluster make negligible contributions to the mea- 
sured DM. The M81 halo DM contribution is estimated° 
to be ~ 15-50 pecm™™ depending on the choice of 
a halo density profile and baryon fraction in the M81 
halo. For the Galactic contribution, models of the Milky 
Way ISM predict DMMY = 30-40 pecm™ in this di- 
rection”? (with a shallow spacial gradient in the re- 
gion, Figure 5). Using these estimates, we limit the 
Milky Way halo contribution to the FRB sight-line, 


which is poorly constrained by current observations” 


2 


to DMNV = DM- DMN —DMicgm -DMmgi $S 32- 
42pecm™*. This is well in line with the models of 
Yamasaki & Totani®, which predict DMN = 10- 


30 pe cm™3, but lower than predicted in some other mod- 
els °°, 

Similar considerations are also applicable for the ro- 
tation measure® RM = -36.9 rad m~? as determined in 
the companion paper by Nimmo et al. (submitted). 
The Galactic contribution’ along the line of sight to 
FRB 20200120E is RM =-17+4radm™~ (Figure 5). 
The contribution of the IGM is most likely minor, leav- 


ing only the MW halo, the M81 halo and the local en- 
vironment of FRB 20200120E as contributing sources. 
The RMs of the MW- and M81-halo are likely small 
(IRMI < 20 rad m7”), constraining RMroca to the range 
[+20, -60] rad m~?. This is comparable to earlier results ae 
for FRB 20180916B (RM = —115radm~*) but three or- 
ders of magnitude lower than the (highly variable) RM of 
FRB 20121102A”. This indicates that models of FRB 
sources do not necessarily require extreme magneto-ionic 
environments, unless the magnetic fields along the line of 
sight to both FRB 20200120E and FRB 20180916B are 
strongly tangled, such that they result in a low net RM. 
More likely though, FRB 20121102A resides in a very dif- 
ferent environment giving rise to its observed properties. 


Data availability 


The datasets generated from the EVN observations and 
analysed in this study are available at the Public EVN Data 
Archive under the experiment codes EK048B, EK048C, 
and EK048F. 


Code availability 


The codes used to analyse the data are avail- 
able at the following sites: AIPS (http: 
//www.aips.nrao .edu/index.shtml), CASA 
(https://casa.nrao.edu), Difmap (ftp://ftp. 
astro.caltech.edu/pub/difmap/difmap.htm1), 


DSPSR (http: //dspsr.sourceforge.net/), 
FETCH (http: //dspsr.sourceforge.net/), 
Heimdall (https: //sourceforge.net/projects/ 


heimdall-astro/), IRAF (http://ast. 


noao.edu/data/software), PRESTO (https: 
//github.com/scottransom/presto), PSRCHIVE 
(http://psrchive.sourceforge.net), and SpS 


(https: //github.com/danielemichilli/SpS). 
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Figure 5: Dispersion measure (DM) and rotation measure (RM) maps around FRB 20200120E: Panel a shows the 
expected Galactic DM contribution (background) according to the YMW16 model disk contribution only, 3° the DM 
of FRB 20200120E (pentagon) and the DMs of known pulsars from the ATNF Pulsar Catalogue in this field circles, !°. 
Panel b shows the physical Galactic Faraday depth Øg background,” , the RM of FRB 20200120E (pentagon) and 
galactic pulsars with a known RM (circles). We assume that the RM° of FRB 20200120E is —36.9 rad m~? (also see 
Nimmo et al. (submitted) for details). 
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Table 4: Notable properties of [PR95] 30244. 


Property Value Reference 
Metallicity log[Z/Zo] =1.74°%0 this work 
Metallicity [Fe/H] “1.83088. this work 
Stellar mass logi[M/Mo] 5 7702 this work 
Effective radius (Re/pc) 3.7 this work 
Age (Gyr) 9.13371, this work 
(u — r) (AB mag) 1.96(2) a 

E(V-B)? 0.201 this work 
c, (kms7!) 22 this work 
Absolute r-band mag. (AB) —8.4 — 

Luminosity distance (Mpc) 3.6 m 


P Milky Way extinction is corrected using a 


reddening map !°°. 


e Using R, = 3.1. 


Table 5: Broadband SDSS filters used to model the SED of [PR95] 30244. 


Effective 
Instrument Filter 
SDSS u 3546 
g 4670 
r 6156 
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Figure 7: MCMC simulation corner plot, where the posterior probability distributions are shown for each of the five 
model parameters along the diagonal panels, and the correlations between model parameter posteriors are shown along 
the columns. Above each probability distribution, the median of the parameter posterior is printed, along with the 
one-o error bars. 
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